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We numerically investigate the phenomenological nonlocal Ginzburg-Landau Hamiltonian for 
two-dimensional superconductors in a strong magnetic field by Langevin equation. We obtain a 
regular vortex lattice which is very near to the square lattice. We calculate the structure factor 
and Abrikosov factor at various temperatures. We also evaluate the specific heat and obtain a cusp 
which indicates the melting of the vortex lattice. 

PACS numbers: 74.20.De, 74.60.Ge 

The study of the high T c superconductor (HTSC) in a magnetic field has renewed the interests in the mixed state. 
The most characteristic behavior of this unconventional superconductor is the large superconducting fluctuation, and 
it is known that the vortex liquid phase appears below the H C 2 line in the H — T phase diagram. The melting 
transition of the vortex lattice has been studied theoretically ever since the discovery of the HTSC. For the 
two-dimensional case, there has been a large amount of studies which investigate the vortex lattice melting transition 
of the Kosterlitz-Thouless type which is accompanied by the dissociation of the dislocation pairs. (|^] Experimentally, 
the vortex lattice melting transition is observed as the first order phase transition by revealing a sharp step of the 
magnetization in both Bi2SriCaCuiO%j r s and Y Ba^Cu^Oi-^. [^,^) The first order melting transition of the vortex 
lattice is already predicted theoretically from the renormalization group analysis based on Ginzburg-Landau model 
for three dimensional system. |Q Although it is not yet certainly clear whether the first order melting transition 
survive in the strong magnetic field phase where the fluctuation effect is very large and the system behaves more like 
two-dimensional. 

On the other hand, the oxide high T c superconductor has a layered structure and it is considered as a quasi- 
two-dimensional system. It is pointed out that the two-dimensional square lattice structure of CuC>2 plane has an 
effect not only on the symmetry of the order parameter but also on the vortex lattice configuration. || Recently 
many experimental evidences which support the notion that the order parameter of HTSC has d x 2_ y 2 symmetry are 
reported. The expansion of the GL model into a form which deals with the d-wave superconductors are thus necessary 
in order to investigate the vortex state in high-T c cuprates. 

In Ref. we have previously investigated a nonlocal GL model starting with following two motivations. One is 
that the nonlocal GL model has a remarkable similarity to the n-component GL model. We find that a parameter of 
the characteristic range of the nonlocal interaction plays a role of n for the n-component GL model. The other one 
is that because of the shortness of the coherence length of HTSC, modification of the quartic interaction of the order 
parameter into a nonlocal form is considered to be reasonable. 
We introduce a phenomenological nonlocal Hamiltonian, 
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where the kernel g(ri,r2) represents a nonlocal interaction between \ip\ 2 - For the symmetric case, we choose 
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where d indicates the range of the nonlocal interaction. In Ref. |9| the high temperature series for the free energy is 
obtained with this Gaussian form nonlocal interaction. 

In this study, we investigate numerically the model Hamiltonian Eq. ([!]) in the case that the kernel g(ri,r2) has 
a four-fold symmetry with Langevin equation. This model is studied recently in Ref. |l0| and they calculate the 
structure factor using the approximation which is valid at moderately low temperatures. We take the same notation 
as Ref. pot , and the Fourier transform of g(r) is written as 



g(k) = exp 
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where k = fc(cos#, sin(9) and n^ 1 = y/l/2eB is the magnetic length. The two parameters C and e determine the 
amount of the effect of nonlocal interaction on the system. We consider the two-dimensional system in a magnetic 
field perpendicular to the xy-plane. We also take the unit h — c = fcg = 1 and the gauge A = (0, xB). 

In a magnetic field, the order parameter ip(r) is expanded with the Landau levels. In order to simplify the numerical 
calculation, we take the lowest Landau level approximation which is valid in a strong magnetic field. Then tp(r) is 
written as 

V(r) = £ a^-VV/^r^exp^y - £{x - q/^ 2 ) 2 } , (4) 
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where q is the momentum in xy-plane, a q is a complex number and L y is the system size. The Hamiltonian is rewritten 
with this V as 
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In our numerical simulation study, the momentum is discretized as 
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where N is the number of vortices in the system with length L x and L y , 

N = \U^ht, (7) 

2tt v ; 
The Langevin equation for the complex order parameter a q = A q + iB q is given by 

dA„ 1 dH 

-df--20A q + " (8) 
dB„ 1 dH 

-i = -23B q ^ (9) 

where 77 is the Gaussian white noise. With these Langevin equations we can calculate a set of a q at each time step, 
and by averaging over the time steps we obtain the statistical average for the values of a q . We denote the discretized 
time step by 5. We rewrite A q = A(I) ,B q = B(I) where I indicates the discretized momenta I = 1, . . . , N. After a 
time evolution of step S, A(I),B(I) becomes A' (I) and B'(I) which is given by 

A'(I)=A(I)-^(2a H A(I) + Si) + Vs v (10) 

B'(I) = B(I)-^(2a H B(I) + S2) + VS V (11) 
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and 
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We set typical time step S = 0.05 and discard first 40000 steps for the equilibration. We put N — 64 vortices in the 
system, and the initial condition is taken as A(I) = B(I) = for all I. We also take the periodic boundary condition 
in the momentum space. 



A(I) = A{mod{I,N)) 
B(I) = B(mod(I,N)) 
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We carry out the numerical simulation for various reduced temperatures y = (T — T C )/AT which is related to an as 

y = a H VV^A (18) 

Fig. 1 shows a configuration of vortices at temperatures y = — 10 and y = —4. The parameters are taken as 
C = 0.01 and e = 0.5 in this calculation. We plot < \tp(x,y)\ 2 > whose averaging is done over 1.6 x 10 5 time steps. 
The white part indicates the large value of < |V>| 2 > and the black spots correspond to the vortices. At y = —10 we 
see a regular rhombic vortex lattice structure instead of the triangular lattice. At a temperature y = —4, although 
the magnitude of < IV'I 2 > is much lowered the rhombic configuration is still maintained. 

The vortex state is probed by the density-density correlation function which is called structure factor. This function 
is defined in the momentum space as 
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which is directly related to a neutron diffraction pattern. In the vortex lattice phase 

A(k) = A <5 k ,G 
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where G is the reciprocal lattice vectors. In the case of triangular lattice, the diffraction pattern shows six-fold 
sharp peaks in the first Brillouin zone. In Fig. 2 we plot A(k) = A(k)/A for three different temperatures. In 
the present nonlocal GL model, the six fold symmetry of the diffraction pattern disappears and a rectangular peaks 
which correspond to the nearby reciprocal lattice points for the rhombic lattice are obtained. The reciprocal lattice 
points including the next nearest ones form a stretched hexagon. This configuration almost does not change over the 
temperature region — 12 > y > —4. The height of the peaks, however, become very low at a temperature y ~ —4 and 
this behavior resembles to that of the diffraction pattern obtained by the local GL model at a temperature just below 
the melting point. In high temperature region y > —4, the peaks are totally disappeared. (Note that this temperature 
is still lower than the T c for the mean field theory.) The ring like diffraction pattern for the vortex liquid phase is not 
obtained in the present model. 

We also calculate the Abrikosov factor 
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where < • • • > denotes the spatial average. In the vortex lattice state (3a is given by 
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where the summation is taken over the reciprocal lattice vectors G. For the triangular lattice (3a = 1-16 and for the 
square lattice Pa — 1.18. At temperature y = —12 we obtain Pa — 1-20 which is larger than the value for the square 
lattice. For the rhombic lattice Pa is evaluated as 



Pa = exp {m 2 cos 2 a + (n — m sin 2 a) 2 } 

^ L cos a 
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where a is the half of the angle which is made by two adjacent sides. For the case of Pa = 1.20, this angle is evaluated 
as a = 45.84° which is very near to the value for the square lattice. Above the temperature y ~ —6, Abrikosov factor 
grows rapidly at about y = —4 and it converges to the value Pa = 2.0 which corresponds to the uncorrelated vortex 
liquid phase or the normal phase. From this behavior we see that in the nonlocal GL model the melting transition of 
the vortex lattice takes place at a temperature near y = —4. 

Next we calculate the specific heat. The specific heat is obtained by taking a second derivative of Eq. (0) with 
respect to the temperature an- 



where 




In Fig. 3 we plot the scaling function for the specific heat C/AC where AC is the specific heat jump at T c in the 
mean field theory. We obtain a cusp in the scaling function of the specific heat at y = —4 at which temperature the 
Abrikosov factor shows a discontinuous jump. This singularity of the specific heat curve has been previously obtained 
by the matrix GL model approach with large- TV limit where the free energy is calculated up to 0(l/./V). In the 
matrix GL model for the large N limit, the specific heat curve shows remarkable agreement with the usual GL model 
both in low temperature and high temperature region. In the current model, however, such agreement is lost except 
the high temperature region. 

In summary, we studied a phenomenological nonlocal GL model numerically using the Langevin equations. We 
obtained the rhombic vortex lattice in the low temperature phase and the corresponding value for the Abrikosov factor 
[3a — 1-20. The cusp in the specific heat suggests that in the nonlocal GL model the superconducting transition reveals 
a totally different feature from that in the original GL model which shows no singurality at the critical point from the 
high temperature expansion series study. [ p"l"[ In order to determine whether the melting transition is discontinuous, 
detailed studies both with numerical and analytical methods are necessary. 

The numerical calculations were made by using the Vector Parallel Processor, Fujitsu VPP500 at RIKEN. 
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FIG. 1. The spatial distribution of < \%/j(x,y)\ 2 > for the system of N = 64 vortices. The black spots correspond to the 
vortices. 



FIG. 2. The diffraction pattern at three diffent temperatrues. We plot only the nearby peaks from the origin k = 0. 

FIG. 3. The scaling function of the specific heat C/AC. The results from the different iteration times are shown for the 
check of equilibration at various temperatures. The cusp at y = —4 is enhanced by taking a longer iteration time. 
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